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ABSTRACT 

We have made a serendipitous discovery of a massive (^ SxlO^^Mo) cD galaxy at z=1.096 in a 
candidate rich cluster in the HUDF area of GOODS-South. This brightest cluster galaxy (BCG) is 
the most distant cD galaxy confirmed to date. Ultra-deep HST/WFCS images reveal an extended 
envelope starting from ~ 10 kpc and reaching ^ 70 kpc in radius along the semi-major axis. The 
spectral energy distributions indicate that both its inner component and outer envelope are composed 
of an old, passively-evolving (specific star formation rate < 10~^Gyr~^) stellar population. The cD 
galaxy lies on the same mass-size relation as the bulk of quiescent galaxies at similar redshifts. The 
cD galaxy has a higher stellar mass surface density (~ M^/R^q) but a similar velocity dispersion 

(~ \/M*/R5o) to those of more-massive, nearby cDs. If the cD galaxy is one of the progenitors of 
today's more massive cDs, its size (R50) and stellar mass have had to increase on average by factors 
of 3.4 ± 1.1 and 3.3 ± 1.3 over the past ~ 8 Gyrs, respectively. Such increases in size and stellar 
mass without being accompanied by significant increases in velocity dispersion are consistent with 
evolutionary scenarios driven by both major and minor dissipationless (dry) mergers. If such cD 
envelopes originate from dry mergers, our discovery of even one example proves that some BCGs 
entered the dry merger phase at epochs earlier than z = 1. Our data match theoretical models which 
predict that the continuance of dry mergers at z < 1 can result in structures similar to those of massive 
cD galaxies seen today. Moreover, our discovery is a surprise given that the extreme depth of the 
HUDF is essential to reveal such an extended cD envelope at z > 1 and, yet, the HUDF covers only a 
minuscule region of sky (r^ 3.1 x 10~*). Adding that cDs are rare, our serendipitous discovery hints 
that such cDs may be more common than expected, perhaps even ubiquitous. Images reaching HUDF 
depths of more area (especially with cluster BCGs at z > 1) are needed to confirm this conjecture. 

Subject headings: galaxies: clusters: general — galaxies: elliptical and lenticular, cD — galaxies: 
evolution 



1. INTRODUCTION 

cD galaxies are the most luminous and most massive 
galaxies in the local Universe and are mostly located 
close to the centers of galaxy clusters. Earlier stud- 
ies have characterized cD galaxies by their excess light 
('envelopes') over the de Vaucouleurs (or r-^/*, de Vau- 
couleurs 1948) profile or Sersic (or r^/", Sersic 1968) pro- 
file at large radii (e.g., Matthews et al. 1964; Schombert 
1988; Graham et al. 1996). More recent studies have 
identified cD galaxies using improved Petrosian (Pet- 
rosian 1976) parameter profiles (Broughet al. 2005; Patel 
et al. 2006; Liu et al. 2008). A large fraction of massive 
nearby brightest cluster galaxies (BCGs) can be classified 
as cD galaxies (Liu et al. 2008). 

The formation mechanism of BCGs/cDs and their 
evolution are not well-understood. To that end, sev- 



eral formation mechanisms including galactic cannibal- 
ism (White 1976; Ostriker & Hausman 1977), tidal strip- 
ping from cluster galaxies (Gallagher & Ostriker 1972; 
Richstone 1976; Merritt 1985) and star formation on 
BCGs by cooling flows (Fabian 1994) have been pro- 
posed. Recent semi- analytic models and simulations of 
galaxy formation suggested that BCGs form in a two- 
phase process: an initial collapse with rapid cooling and 
star formation at high redshift is followed by later (z < 1) 
growth through multiple dissipationless (dry) mergers of 
pre-existing progenitors (e.g., Dc Lucia & Blaizot 2007; 
Ruszkowski & Springcl 2009; Naab et al. 2009; Laporte 
et al. 2012). Very recently, Laporte et al. (2013) showed 
that BCGs can form through dissipationless mergers of 
quiescent massive z = 2 galaxies in their simulations in 
a ACDM universe, which likely indicates that BCGs en- 
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tered the dry merger phase at epochs ear Her than z = 1. 
Interestingly, they found eD signatures among their sim- 
ulated BCGs by z ~ 1 (Laporte et al. 2013). Recent 
observational studies also provide evidence that BCGs 
build up a large part of their stellar masses via dry merg- 
ers at z < 1 (e.g., Liu et al. 2008,2009; Tovmassian & 
Andernach 2012; Lidman et al. 2012). However, other 
studies argued that BCGs evolve minimally since z ~ 1 
(Whiley et al. 2008; Stott et al. 2011) and mergers 
are unlikely to be important for the growth of BCGs at 
z < 1 (e.g., Ascaso et al. 2011). Our early work (Liu et 
al. 2008) showed that BCGs with dry mergers evolve to 
possess light profiles similar to those of local cD galaxies 
by depositing stellar plumes into the BCG outskirts. 

Although the origin of cD galaxies, and the contribu- 
tion of their envelopes to the intra-cluster light, are still 
not completely clear (e.g., Gonzalez et al. 2005; Zibetti 
et al. 2005), the extended stellar haloes of cD galax- 
ies to surface brightness fi{r) < 27 magarcsec"^ in the 
rest-frame Sloan r-band are likely to come from galaxies 
themselves (Tal & van Dokkum 2011): the intra-cluster 
light has much lower surface brightness and only domi- 
nates at large radius (r > 80—100 kpc; for detailed dis- 
cussions, see Zibetti et al. 2005; Bernardi et al. 2007; 
Lauer et al. 2007; Tal & van Dokkum 2011). Even 
so, it is difficult to identify faint cD signatures in the 
high-redshift (z > 1) galaxies due to the cosmological 
surface brightness dimming. Therefore, it is still un- 
known whether cD galaxies can form at high redshift. 
Best et al. (1998) claimed a detection of cD-like signa- 
tures in two high-redshift SCR radio galaxies (3C 22 at 
z=0.936 and 3C 41 at z=0.795) based on shallow (one or- 
bit) HST/WFPC2 imaging data. However, their results 
were not confirmed by deeper imaging observations. 

In this paper, we report the serendipitous discov- 
ery of a massive cD galaxy (GOODS-S J033237.19- 
274608.1, named after GOODS-S and its J2000 coordi- 
nate: RA=53.1549542, DEC=-27.7689028) at z=1.096 
(see §5) in a candidate rich cluster in the HUDF area of 
GOODS-South. We use its ultra-deep HST/WFC3 mo- 
saic image to ~ 30magarcsec~^ to confirm the existence 
of an extended cD envelope characterized by distinct Pet- 
rosian parameter profiles. We also investigate the stellar 
population and environment of this distant cD galaxy, 
and explore its possible evolution to today by compar- 
ing its structural parameters to those of more-massive, 
nearby cD galaxies. Throughout the paper we adopt a 
cosmology with a matter density parameter fim = 0.3, a 
cosmological constant 51 a = 0.7 and a Hubble constant 
of Hq = 70kms~"'^Mpc^^. All magnitudes are in the AB 
system unless otherwise stated. 

2. IMAGE DATA 

The z - 1.1 cD galaxy, GOODS-S J033237.19- 
274608.1, is located in the Hubble Uhra Deep Field 
(HUDF, Beckwith et al. 2006) region of the south 
field of the Great Observatories Origin Deep Survey 
(GOODS-S, Giavalisco et al. 2004). The HUDF is 
covered by the HST/WFC3 F105W(y), F125W(J), 
and F160W(ii^) bands imaging from G. Illingworth's 
HUDF09 program (GO 11563, Bouwens et al. 2010) 
as well as the Cosmic Assembly Near-IR Deep Ex- 
tragalactic Legacy Survey (CANDELS, Grogin et al. 
2011; Koekemoer et al. 2011). The HUDF has also 



been observed by the HST/ACS F435W(B), F606W(y), 
F775W(i), and F850LP(z) bands imaging from Beck- 
with et al. (2006) and GOODS (Giavalisco et al. 2004). 
We made use of all HST observations published prior 
to June 2012 in the HUDF to maximize our ability 
to detect the low surface brightness envelope of this 
cD galaxy at large radius. These HST mosaics reach 
the 5tT limiting depths (O'.'35-diameter aperture) of 29.7, 
30.1, 29.9, 29.4, 29.6, 29.9, and 29.9 AB in B, V, i, 
z, F, J, and 7J-bands. Low- resolution UV (U-band 
from both CTIO/MOSAIC and VLT/VIMOS) and IR 
(VLT/ISAAC Ks and Spitzer/IRAC 3.6, 4.5, 5.8, 8.0 
^m) images are also available for this galaxy. Uniform 
multi-wavelength photometry from the UV to 8.0 fxm 
with the resolution and blending of each band taken into 
account is provided by Guo et al. (2013). 

3. IDENTIFICATION OF CD SIGNATURE 

To quantitatively determine how a surface brightness 
profile deviates from a de Vaucouleurs or Sersic profile, 
we use an objective method based on the Petrosian rj(r) 
profiles (Liu et al. 2008), defined as 



■>l{r) = m(0- (m(0>: 



(1) 



where /x(r) is the surface brightness in magnitudes at ra- 
dius r and (fJ-ir)) is the mean surface brightness within r 
(Petrosian 1976). We also calculate the gradient of r]{r), 
7(r) = dri{r)/d\og{r). Liu et al. (2008) have shown that 
there is always a distinct signature of a plateau in the 
Petrosian r]{r) profiles and a valley in the "f(r) profiles 
for cD-like BCGs with an extended stellar halo (see their 
Fig. 5). Such signatures are not present in normal ellip- 
tical galaxies, which are usually well-fit by a r^/" profile. 
We label the radius of the minimum in the j{r) profile as 
Rmin- We also define a depth, Z^min, which is the difi'er- 
ence between the minimum and maximum of 7(r) outside 
twice the seeing radius and Rmin- Rmin is, by definition, 
the transitional radius (TR) between the inner compo- 
nent and outer halo in a cD-like galaxy. Our method 
indicates that the surface brightness deviations of over- 
all BCGs from the single Sersic profile does not have any 
sharp transitions. However, the more significant the de- 
viations of the surface brightness profiles from a single 
Sersic profile (i.e., the more prominent the stellar halo in 
the outskirts), the more distinct the signatures (plateau 
and valley) become statistically (see Liu et al. 2008 for 
details) . 

Fig. 1 is the pseudo-color composite image of GOODS- 
S J033237. 19-274608.1 from the ACS F-band, WFC3 J- 
and i/-band images. The transitional isophote and the 
maximum measured isophote are marked with red el- 
lipses. Panels (a)-(c) of the Fig. 2 show the observed 
surface brightness profile, the corresponding Petrosian 
ri[r) profile and 7(r) profile in the J-band, all of which 
are measured from the mosaic image of CANDELS/Deep 
+ HUDF09 (5 orbits + 34 orbits). AU radial profiles 
use the equivalent radius of an ellipse, rab=vab, where 
a and b are the lengths of semi-major and semi-minor 
axes of the ellipse. Contaminating objects detected with 
counts above 1 sigma noise have been masked during our 
measurement following Liu et al. (2008). The F125W 
(J) band is used because it is closest to the rest-frame 
Sloan r-band at this redshift. It is clear that the surface 




Fig. 1.— The pseudo-color image (0.5 Mpc X 0.5 Mpc) of GOODS-S J033237. 19-274608.1. North is up and east is to the left, 
transitional isophote (inner) and the maximum measured isophote (outer) are marked with red ellipses. 



The 



brightness profile of GOODS-S J033237. 19-274608.1 de- 
viates significantly from a single de Vaucouleurs (r^'^) 
profile or Sersic (r^/", n ^ 3.2 here) profile, and the 
r]{r) profile shows a plateau and the j{r) profile shows 
a valley around the transitional radius. The transitional 
radius i?min is 9.7 kpc and the corresponding -Dmin value 
is 1.35. 

We investigate the position of GOODS-S J033237.19- 
274608.1 on a plot of Anin versus i?min compared to 
a representative sample of nearby BCGs in Liu et al. 
(2008), as shown in the panel (d) of Fig. 2. Liu et al. 
(2008) statistically showed that the larger Dmin and the 
transitional radius Rmin of a BCG are, the more sig- 
nificantly extended its envelope appears to be. BCGs 
with Dniin larger than 0.8 are more likely to be classi- 
fied as cD galaxies. Note that -Dmin — 0-8 is the me- 
dian value for all nearby BCGs. Our galaxy, GOODS-S 
J033237. 19-274608.1 at z = 1.096, has £)„!„ = 1-35 and 
is therefore termed a cD galaxy. The stellar envelope 
of this high-redshift cD galaxy starts to dominate from 
Tab = 9.7 kpc (the corresponding semi-major axis radius 
is a = 11.6 kpc) and extends to 65.3 kpc (a = 70.12 kpc). 

4. THE SIZE AND STELLAR MASS 

It is unreasonable to characterize the sizes of cD galax- 
ies by half-light radii (effective radii) from single Sersic 
(n = 4 for de Vaucouleurs) fits since the surface bright- 
ness profiles of cD galaxies usually show strong devi- 
ations from single Sersic profiles. We opt to use the 
real isophotal half-light radius to characterize the size 
of this cD galaxy following Liu et al. (2008). The de- 



rived half-light radius of GOODS-S J033237. 19-274608.1 
within the maximum measured isophote in the J-band is 
only ~ 5.25 kpc. In order to estimate the effects of point- 
spread-function (PSF) on the size measurement, we first 
used the GALFIT package (Peng et al. 2002) to fit a 
conventional two-dimensional model, convolved with the 
real PSF, to the observed image. A double-Sersic model 
was applied to determine the best-fit model image. We 
then added the residual image to the best-fit model image 
(PSF deconvolved) to correct the uncertainty of model 
fitting. It should be noted that the residual image is not 
deconvolved for PSF. However, this does not significantly 
affect the accuracy of our method since the flux ratio of 
residual-to-observed images is only ^ 3%. We did the 
isophotal photometry and half-light radius measurement 
for the 'residual-corrected' model image in the same way 
as for the observed image. The obtained half-light ra- 
dius is ^ 5.19 kpc, which shows the PSF effect is not 
significant. 

We followed Guo et al. (2012) to fit spectral energy 
distributions (SEDs) composed of homogeneous aper- 
ture magnitudes from the UV to 8fim to estimate the 
stellar mass. We used an updated version (CB07) of 
Bruzual & Chariot (2003) stellar population synthesis 
models with a range of stellar population properties and 
dust extinctions. We assumed models with a solar nietal- 
licity, a Salpcter (Salpeter 1955) initial mass function 
(IMF), and a Calzetti et al. (2000) law. The derived 
stellar mass is lO^'^^io.os^^^ j|- becomes -- 5x10"Mq 
j-2Qii.7i±o.08]y[^-j a,fter accounting for the small difference 
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Fig. 2. — The panels (a)-(c) show the surface brightness profile (a), the corresponding Petrosian ri{r) profile (b) and 7(r) profile (c) in 
the J-band, respectively. The vertical dotted lines mark twice the PSF radius. The vertical solid lines mark the transitional radius (TR). 
The solid curves in panel (a) are the best-fit de Vaucouleurs (black) and Sersic (magenta) profiles fit to the data between twice the PSF 
radius and maximum measured radius. Panel (d) illustrates the position of GOODS-S J033237. 19-274608.1 on the relation of -Dmin versus 
Rmin for nearby BCGs in Liu et al. (2008). Nearby BCGs are shown with circles and GOODS-S J033237. 19-274608.1 is marked with a 
red star. Symbol sizes are proportional to the stellar masses that are derived from the SED fitting (see §7). The horizontal dashed line 
shows the median value (0.8) of D„iin in nearby BCGs 



between the used aperture magnitude and isophotal mag- 
nitude within the maximum measured elhpse in the J- 
band. 

5. OPTICAL SPECTROSCOPY AND STELLAR 
POPULATION 

We show the optical spectroscopy of GOODS-S 
J033237.19-274608.1 in the top panel of Fig. 3, which 
is observed by the VLT/F0RS2 instrument in the 
ESO/GOODS spectroscopic program. The grating used 
has a scale of roughly 3.2 A/pixel and a nominal res- 
olution of A/AA = 860. The whole spectral range is 
6000- lOSOOA (see Vanzella et al. 2005, for details). The 
derived spectroscopic redshift of GOODS-S J033237.19- 
274608.1 is 1.096±0.001. The spectroscopy of this galaxy 



shows weak emission lines (e.g., [OII]3727) and a weak 
H(5 absorption line, indicating that it is dominated by an 
old stellar population. 

We fitted the SEDs of the inner component (r < 
Rmin) and outer envelope {Rmm < r < 30 kpc) to 
investigate each of their stellar populations. The un- 
masked imaging data in seven high-resolution HST bands 
{B,V,i,z,Y,J,H) were used. The cut off at 30 kpc for the 
outer region is to ensure high signal-to-noise ratios for 
all bands at large radius. The results are shown in the 
bottom panels of Fig. 3. The derived specific star for- 
mation rates (sSFRs) for the inner component and outer 
envelope are 3.7 x lO^^Gyr"^ and 1.4 x 10~'*Gyr~^, re- 
spectively. The sSFRs, i.e. the SFR per unit stellar 
mass, were calculated by dividing the SFRs by the stellar 
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Fig. 3. — The top panel is the optical spectrum of GOODS- 
S J033237.19-274608.1 observed by the VLT/FORS2 instrument. 
The bottom panels are the best-fit SEDs composed of high- 
resolution HST bands for the inner component (bottom left) and 
the outer envelope (bottom right), respectively. 



masses, which were derived from the SED fittings sepa- 
rately. The sSFRs of both components are consistent 
with being passively-evolving stellar populations. 

6. THE ENVIRONMENT 

It is well known that cD galaxies in the local Uni- 
verse are usually BCGs and are mostly found in the cen- 
ters of rich clusters. We now investigate whether this 
high-redshift cD galaxy is also located in a cluster en- 
vironment. Salimbeni et al. (2009) have found three 
group/cluster structures around the photometric redshift 
Zp ~ 1.05 in the entire GOODS-S field. The projected 
density isosurfaces of these three structures are shown in 
the top image of Fig. 4 (also see panel c in Salimbeni et 
al. 2009). It is obvious that the discovered z ~ 1.1 cD 
galaxy is very close to the projected density peak of the 
structure with ID=7. Salimbeni et al. (2009) estimated 
that this structure has a photometric redshift Zp ^ 1.04, 
1200 ~ 1-2 — l.SMpc and 60 member galaxies. Our cD 
galaxy is the most massive galaxy within the redshift gap 
of Azp = 1.04 ± 0.1 in this structure. We searched for 
other galaxies with secure spectroscopic redshifts (qual- 
ity flag = 1 or 2) in the literature (Vanzella et al. 2005, 
2006; Balestra et al. 2010; Cooper et al. 2012) within a 
projected distance of 1 Mpc from the density peak of this 
structure. The distribution of redshifts (radial velocities) 
of these galaxies (with a bin size of 0.001, corresponding 
to typical errors of spectroscopic redshifts) near the es- 
timated structure redshift is shown in the bottom panel 
of Fig. 4. We iteratively fitted this distribution with a 
Gaussian function and rejected all outliers which lie be- 
yond 3 sigma of the best-fit Gaussian distribution. The 
final best fit includes 21 objects all within 3 sigma of 
the best-fit Gaussian, with an inferred velocity disper- 
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Fig. 4. — The top image illustrates the position of GOODS- 
S J033237. 19-274608.1 (red cross) on the WFC3 J-band image of 
the entire GOODS-S field. The projected density isosurfaces for 
structures at Zp ~ 1.05 identified by Salimbeni et al. (2009) are su- 
perimposed on it. Yellow crosses indicate the density peak of each 
structure, and the number of each is the ID of the structure given 
by Salimbeni et al. (2009). The big yellow circle marks the pro- 
jected range within 1 Mpc from the density peak of the structure 
with ID=7. The bottom panel shows the distribution of redshifts 
(radial velocities) of galaxies with secure spectroscopic redshifts 
near the cluster redshift and within the marked range (yellow cir- 
cle). The curve is the final best- fit Gaussian distribution. Red 
histograms represent the data used during the final fitting. Black 
histograms represent galaxies which are likely outliers rejected by 
our iterative fittings. The inset is designed to show the distribution 
of spectroscopically confirmed members more clearly. 



sion of 810±120km/s. We derived a precise redshift lo- 
cation of z = 1.0965 ± 0.0027 for this structure. These 
values are consistent with the results (797±138km/s and 
z = 1.0968 ± 0.0027) estimated using the biweight esti- 
mators, computed with the Rostat package (Beers et al. 
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1990). The redshift of the cD galaxy places at the cen- 
ter of the cluster redshift distribution. These properties 
collectively suggest that GOODS-S J033237.19-274608.1 
is both likely to be in a cluster environment and is also 
located close to the cluster center. Future spectroscopic 
observations can confirm more cluster galaxies and quan- 
tify the richness of this potential rich cluster. 

7. SUMMARY AND DISCUSSION 

We have made a serendipitous discovery that the most 
massive (~ SxlO^Mg) galaxy (GOODS-S J033237.19- 
274608.1), i.e., the BCG, in a candidate rich cluster at 
z=1.096 in the HUDF area of GOODS-S is a cD, with an 
outer envelope extending in radius to ~' 70 kpc along the 
semi-major axis. This is the most distant, massive cD 
galaxy confirmed to date. The HUDF image is unique in 
being the deepest and, thus, as seen in Fig. 5, such a faint, 
extended, stellar envelope would not have been otherwise 
visible. Moreover, since the HUDF only covers less than 
5 square arcminutes (i.e., 1 part in 30 millionth of the 
sky) and BCGs are expected to be very rare at z > 1, our 
discovery of even one such BCG being a cD was a total 
surprise and is certainly fortuitous by any expectation. 
Beyond the following discussion of the implications for 
the early and later formation of some cDs, our discovery 
hints that many other z > 1 BCGs might be found with 
such extended cD halos, if more areas or clusters were 
covered with deep enough images. 

A clue to the early formation of cDs comes from the 
SEDs of both the inner component and outer envelope of 
this cD galaxy; the SEDs match those of old, passively- 
evolving stellar populations. Assuming that such cD en- 
velopes originate from dry mergers (see Liu et al. 2008), 
our discovery of even one such example at z > 1 proves 
that some BCGs entered the dry merger phase at epochs 
earlier than z = 1, a result consistent with recent simu- 
lations by Laporte et al. (2013). 

Next, we explore the possible evolution of our z > \ 
cD galaxy to today by comparing its structural parame- 
ters (i.e., mass, size, and velocity dispersion) to those of 
the more massive cDs in the local Universe. The stellar 
masses of local cDs were derived using the same CB07 
models with the same Salpeter IMF, solar metallicity, 
and a Calzetti extinction law (Calzetti et al. 2000) as 
adopted for deriving the mass of our z > 1 cD galaxy. 
The SEDs used Petrosian magnitudes in ugriz bands 
from the Sloan Digital Sky Survey and the extrapolated 
magnitudes in JHKg bands from the Two Micron All 
Sky Survey Extended Source Catalog (2MASS XSC). 
The final masses were re-scaled by adding the differences 
between the Sloan r-band Petrosian magnitudes and 
isophotal magnitudes measured within 25magarcscc~^. 

The r-band isophotal half-light radius of nearby BCGs 
within 25 magarcsec"^ were measured by Liu et al. 
(2008). To make the comparison at the same rest- 
frame bandpass, we applied the best-fit global SEDs of 
GOODS-S J033237. 19-274608.1 to convert the observed 
J-band surface brightness profile to the rest-frame Sloan 
r-band. The surface brightness profile has also been cor- 
rected for the cosmological dimming of 10 log (1+z). 
The corresponding stellar mass and half-light radius of 
GOODS-S J033237.19-274608.1 within 25magarcsec"2 
in the rest-frame r-band are 10"-^^=^°°^Mq and 5.02 ± 
0.12 kpc, respectively. These errors include both random 



errors and errors from the band shifting and PSF effects. 
There are 64 nearby BCGs with log M, > 11.67 out 
of the 85 BCGs in Liu et al. (2008), 33 of which have 
^min > 0.8 and can be classified as cD galaxies with high 
probability. 19 of those 33 nearby cDs have the measured 
velocity dispersions available in Liu et al. (2008). We de- 
rived a velocity dispersion of 333±83km/s for GOODS-S 
J033237. 19-274608.1 from its optical spectroscopy, which 
is consistent with the value (324 ± 32km/s) estimated by 
van der Wei et al. (2005, corresponding to ID CDFS-18 
in their paper). We show the stellar mass surface den- 
sities ('-^M^/Rjq)^, sizes (R50), M*/R5o (sometimes used 
to infer velocity dispersion, see Franx et al. 2008) ^, and 
velocity dispersion (cr) as a function of stellar masses of 
this high-rcdshift cD galaxy and more-massive, nearby 
cDs in Fig. 6, respectively. 

The distant cD galaxy clearly has a higher stellar mass 
surface density (i.e., more compact) but a similar veloc- 
ity dispersion (~ ^M,/R5o) to those of more-massive, 
nearby cDs. It lies, however, on the same mass-size rela- 
tion as the bulk of quiescent galaxies at similar redshift 
(1.0 < z < 1.5), as derived by Newman et al. (2012). 
The more-massive, nearby cDs, in contrast, lie on the 
z ^ mass-size relation (within the scatter) derived by 
Shen et al. (2003) who used Petrosian-bascd quantities. 
To evolve the distant cD galaxy into one of the more 
massive nearby cDs, its size and stellar mass have to in- 
crease, on average, by factors of 3.4 ± 1.1 and 3.3 ± 1.3 
over the past ~ 8 Gyrs. 

As mentioned in the introduction, recent studies from 
both numerical simulations and observations show that 
BCGs accumulate a large part of their stellar mass via 
dry mergers at z < 1. Ruszkowski & Springel (2009) 
first investigated the effect of dry mergers on the scaling 
relations of BCGs in simulations. Laporte et al. (2012) 
argued, however, that the stellar masses and sizes of the 
final merger remnants derived from the Ruszkowski & 
Springel simulations were too massive (~ lO^'^Mo) and 
large {^ 100 kpc) compared to real BCGs at z ^ 0. They 
thus revised the simulations of Ruszkowski & Springel by 
including the effects of initial compactness on the subse- 
quent evolution of the mass-size relation of the BCGs 
(Laporte et al. 2012). The revised simulations showed 
that the size evolution of the most extended BCGs (e.g., 
cDs) at high redshift scale linearly with stellar mass 
(i?5oocM*). Nipoti et al. (2009) also investigated the dry- 
merging evolution scenario based on a set of simulations 
and showed that major and minor dry mergers increase 
half-light radius with stellar mass as i?5oocM^°^^°-^^ but 
have little effects on velocity dispersion (crocM^^^^"-^^). 
The reasonable agreement between our observations and 
theoretical predications mentioned above shows that the 
continuance of dry mergers at z < 1 can result in struc- 
tures similar to those of massive cD galaxies seen today 
and thus support the view that dry mergers play an im- 
portant role in the late evolution of cD galaxies. 

^ The stellar mass surface density is actually is M*/27rRg(j, but 
we omit the constants. 

^ The true stellar velocity dispersion is a^ oc GM/R50, where M 
is the total mass including stars, gas, and dark matter. Pranx et 
al. (2008) provide a value of the coefficient through the fitting of 
a sample of SDSS galaxies: a^ = 0.3GM,/R5o. 
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Fig. 5. — Comparisons of the J-band images (top) and surface brightness profiles (bottom) between the CANDELS/Deep only and the 
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